Coherent Stimulated X-ray Raman Spectroscopy; Attosecond Extension of RIXS 
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Spontaneous and stimulated resonant inelastic X-ray Raman scattering signals are calculated 
using the Keldysh-Schwinger closed-time path loop and expressed as overlaps of Doorway and Win- 
dow electron-hole wavepackets. These are recast in terms of the one-particle Green's functions and 
expansion coefficients of configuration interaction singles for valence excitations, which can be ob- 
tained from standard electronic structure codes. Calculation for many-body states of ground and 
core-excited system is avoided. 



Resonant nonlinear spectroscopy in the X-ray regime 
may become possible by new bright ultra-fast sources 
U B S BL HI • The theoretical formulation of nonlinear 
spectroscopy with attosecond X-ray pulses is of consider- 
able interest [f| 0, B • The picosecond optical pump / x- 
ray probe technique has been used to study photophysical 
and photochemical molecular processes [T(| [U, HH, [l3[ . 
An all- X-ray pump-probe experiment with attosecond X- 
ray pulses has been proposed in Ref. @ . The pump pulse 
interacts with the system to create a valence excited state 
wave-packet which evolves for a controlled delay time r 
when a second probe pulse interacts with it. r is not lim- 
ited by the core-hole life-time. The dependence of this 
coherent stimulated X-ray Raman signal (CXRS) on the 
delay time r carries information about valence excited 
state dynamics. The pump-probe signal may be recast 
in the Doorway /Window representation of optical non- 
linear spectroscopy [3| . In Ref. Q these were computed 
within the equivalent core single Slater determinant ap- 
proximation. Here we express them using the single-body 
Green's functions. Thus avoiding the explicit computa- 
tion of the many electron core excitations. 

Green's functions have been extensively used to study 
X-ray absorbtion fine structure (XAFS) [H [H, EE El, 
[TqI . [20| . The formalism is well developed and incorpo- 
rates intrinsic and extrinsic losses fl9| coming from the 
many-body interactions (electron-phonon, electron-hole 
pairs etc.). These provide a high level, yet practical, 
approach that goes beyond the density-functional the- 
ory considered in Ref. [8]. By invoking the "sudden 
approximation" whereby a core-hole is created and de- 
stroyed instantaneously we express the signal in term of 
one-particle Green's functions which depend on the core- 
hole parametricalrypjj. 

We compare CXRS with resonance inelastic X-ray Ra- 
man scattering (RIXS) which is a common frequency- 
domain technique used in the study of core excited states 
in solids and molecules [H [H, HI H- The RIXS signal 
can be described by the Kramers-Heisenberg expression 
[23L [25| as is done for valence excitations in the optical 
regime [l4| . 
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here io\ and 0J2 are the incoming and scattered photon fre- 
quencies, a and c denote the valence iV-electron ground 
state and singly excited states and e is the excited state 
with one core-hole and N + 1 valence electrons. B ea and 
B ce are matrix elements of dipole operator. 

Displaying this signal vs. u>\ — L02 reveals the valence 
transitions (cj oa ). The core- hole resonances uj ea in the 
transition amplitude A ca are typically much broader and 
less resolved due to the large core-hole lifetime contribu- 
tion to r ea is ~ 0.05 eV [25j]. CXRS is a closely related 
technique, performed with a pair of attosecond pulses. 
We show how that how the signal can be recast in a form 
resembling Eq.([T]), making it a natural time-domain ex- 
tension of RIXS. Varying the envelopes of both pulses 
and their delay offer a much higher degree of control of 
the signal. Both signals can be described as a valence 
electron-hole wavepacket (a doorway state) prepared by 
the pump beam with high spatial and temporal resolu- 
tion. This wavepacket is probed by projecting it into a 
second (window) wavepacket prepared by the probe [2(| . 
This wavepacket can be visualized in real space when 
the valence excitations are treated at the configuration- 
interaction-singles (CIS) level. 

Equation ([I]) can be derived in the time domain us- 
ing the closed-time-path loop (CTPL) diagram shown in 
Fig. 1. This diagram may then be modified to repre- 
sent CXRS. In this diagram time runs clockwise and all 
the interactions are ordered along the loop. In physical 
time, however, only interactions on the same branch of 
the loop are time-ordered with respect to each other; In- 
teractions on different branches are not time-ordered|27j . 
The correlation- function expression for the RIXS signal 
is readily obtained from Fig. 1 by assigning each inter- 
action an excitation operator (B or B') according to the 
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FIG. 1: CTPL diagram representing the RIXS [Eq. ©]. 
Time runs on the loop clockwise starting from bottom of 
the left strand. Interactions on each branch are time-ordered 
with respect to each other. Interactions in different branch 
are not time-ordered. The incoming and detected modes are 
denoted by the indices 1 and 2, respectively, si, S2 and S3 
are positive time intervals along the loop. By ordering the 
interactions in different branches, Fig. 1 can be decomposed 
into three diagrams corresponding to three different Liouville 
space pathways [14{. However, this will not be necessary here. 

rules given in Ref. 27j. We then get 

/oo />oo />oo 

ds 2 / dsi / ds 3 
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(B m {t - s 3 )Bl(t)B n {t ~ s 2 )B\{t -sx- s 2 )). (3) 

The operator B^(B m ) creates (annihilates) a core- 
hole/ valence-electron pair. 

B ,ri ^ ^ Mira^i^m) ^ ^ UrnjCmCj (4) 

i i 

Here c m annihilates electron at the mth core orbital and 
cj creates an electron at the ith valence orbital. These 
are Fermi operators. lo\ and lo 2 are the carrier frequen- 
cies of the two pulses and fjbi m is the dipole matrix ele- 
ment between the ith valence orbital and nth core orbital. 
Ej(r,t),j — 1,2 are complex field envelopes 

Ej(r, t) = Eje ik ^ r - iuj ^ + £* e -^- r +^* (5) 

The time dependence of exciton operators in Eq. ([3]) 
is given by the free molecular Hamiltonian (no field), 
B m (t) = e tHt B m e~ lHt . By inserting the identity opera- 
tor expanded in the many-body states of the 
molecule, between the exciton operators in the middle of 
Eq. ©and assuming stationary field Ej(t) = lwe can 
carry out the s 2 integration and we immediately recover 
the Kramers-Heisenberg (KH) expression ([l) for the sig- 
nal. 

The single-particle many-body Green's functions pro- 
vide a standard tool for computing X-ray absorbtion near 
edge spectra (XANES)[19j. Computer codes based on the 



GW approximation developed for these Green's functions 
are broadly applied to molecules and crystals [l5l IHI [l7l. 
[l8j . We next show that the signal can be approximately 
expressed in terms of these Green's functions. The cor- 
relation function in Eq. ([3]) can be recast as, 

(B m (t - s 3 )Bl(t)B n {t - s 2 )Bt(t - si - s 2 )} 
= (a\B m U(-s 3 )BlU(s 2 )B n U(s 1 )Bi\a) (6) 

where U (s) = e lHs is the time evolution operator and H 
is free molecular Hamiltonian. We set the ground state 
energy u> a = 0. 

We next define the following projection-operator in the 
N + 1-valence electron/ 1-core- hole space. 

V m = ^2c\c m \a){a\c\ n Ci (7) 

i 

where \a) is the ground many-body state with N- valence 
electrons. V m selects a sub-space of the full iV+l-electron 
space which includes single valence electron-hole pair ex- 
citations. By inserting the projection operator ([7]) twice 
inside the average on the r.h.s. of Eq. ©, we obtain 
an approximate expression for the correlation function. 
Using Eqs. ^ and ([7]), it factorizes into a product of 
three correlation functions. 

(B m (t - s 3 )Bl(t)B n (t - s 2 )Bl(t - si - s 2 )} 
« (a\B m U(-s 3 )V m BlU{s 2 )B n V n U(si)Bt\a) 

= ^imA i mfeA i jnA i ni( a l c L C ^(~ S 3)4 C ™l a ) 
ijkl i'j< 

x {a\ciic\U{s2)c ] c^\a){a\c\ l Cj>U{si)c\c n \a). (8) 
Since the core-holes are highly localized on the parent 
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FIG. 2: Three level sequential transition dipole scheme used 
for pump-probe signal, a and c are valence ground and excited 
states, and e is the valence (JV + l)-electron excited state in 
presence of one core-hole. 

atom, their dynamics is very slow compared to the va- 
lence electrons and may be ignored. We shall ignore this 
dynamics and treat the core-hole indices as fixed param- 
eters. 

We next introduce the following set of many-body 
states 

\ X (i,j))=c i c]\a). (9) 
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These represent one electron-hole pair excitation state 
of the valence iV-electron system with one electron-hole 
pair. 

We define the one-electron Green's function computed 
in the presence of a core-hole at m. 



G^\t,t')=-i{Tc i {t)c]{t')) 



t/V\ 



(10) 



where T is the time-ordering operator which rearranges 
a product of operators in increasing order in time from 
the right to the left and (-) m represents a trace over 
iV-electron ground state of the valence in presence of 
Coulomb potential due to a core-hole at m. 
We also define the one-sided Fourier transform 



G$ n) ("i) = 



dt e- iult G${t) 



(11) 



where G R is the retarded Green's functions [28l|. 

Using Eqs. © and (flU)) , the correlation-function © 
can be recast as 

(B m (t - s s )Bl(t)B n (t - s 2 )Bl{t - Sl - a a )> = 

ijkl i'j' 

47 )t (^)G^ ) ( Sl )(x( l ',fc)|e-^|x(j,/)). (12) 

The many-body eigenstates |c) (with energy e c = u> ac ) 
of the molecular Hamiltonian H at the configuration in- 
teraction single (CIS) level are given by the linear com- 
binations of single electron-hole pair states. 



(13) 



where / c; ,j are expansion coefficients. \c) constitute an 
orthonormal single electron-hole pair basis set, 
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Inserting this identity into the last term in Eq. fl2|) . and 
substituting it in Eq. ^ gives 

Srixs (^11^2) — 

^ ^ ^ fc^'kfc-J'jf-^imf-^rnkf-^jnf^nl 
ijkl i'j' mnc 
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where uj ac = uj a — uj c . Equation (|15p correspods to the 
RIXS signal obtained earlier using the Keldysh Green's 
functions approach [Eq. 47 in Ref. [30]]. However 
unlike in Ref. [3(| , the RIXS signal derived here takes in 



to account of the shape of incoming pulse which provides 
a better control on the signal. 

Assuming stationary field envelopes, the S2 integral in 
Eq. (fT5|) reduces to a Dirac-delta function and Eq. (fT5j) 
takes the KH form (TTJ with the transition amplitude 

A ca {u)t) = ^ fc;]j'Gf,\ n) (U^fljn^rd (16) 
nljj' 

Equation (frBJmay be interpreted as follows: The N- 
electron ground state \a) is excited by the incoming X- 
ray beam, creating a core excited state with (N + 1) 
valence electrons. This state evolves during the short 
time window permitted by the core-hole lifetime. This 
evolution in the presence of a core-hole is described by 
the frequency-dependent Green's function. Finally this 
excited state is transformed by a dipole transition to a 
singly excited valence ground state \c). 



I. X-RAY PUMP-PROBE SIMULATED RAMAN 
SIGNAL 

We assume high probe intensity so that stimulated 
emission is dominant and spontaneous emission can be 
neglected. The signal is given by the difference be- 
tween the transmitted intensities of the probe (£2) with 
and without the pump (fci). Four CTPL diagrams con- 
tribute to the signal for our two-band model (Fig. 2). 
When the pump and probe pulse envelopes arc well sep- 
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FIG. 3: Two CTPL diagrams that contribute to the pump- 
probe signal for a two valence level scheme sketched in Fig. 
2. 

arated in time, the process may be separated into three 
steps. First, two interactions with the pump create a 
wavepacket of electron- hole pair states |c) in the valence 
iV-electron system. This wave-packet evolves during the 
delay r and modulates the absorbtion of the probe. The 
two loop-diagrams shown in Figs. 3a and 3b contribute to 
the signal. These diagrams represent two different cor- 
relation functions of the exciton operators. Details are 
given in Appendix [XJ The signal can be expressed in the 
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Doorway/ Window form 

scxrs(t) = 



[DacWc 



+ C.C 



The pump prepares a Doorway wavepacket (|A6 



(17) 



|D) = 5> oc |c) (18) 

C 

and the probe creates a Window wavepacket (|A7j) 



\W) = Y,W ac \c). 



(19) 
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FIG. 4: Schematic representation of Doorway and Window in 
Eqs. (|21]| and (ggj) . 



The signal is given by the overlap of these wavepackets. 
Note that the signal in (fTT|) is given for r > (pump 
interacts before the probe pulse). However, in order to 
express the CXRS signal in a form similar to the RIXS 
(|2ip. for r < 0, we define the CXRS signal to be the same 
as in Eq. fTTJ). 

Fourier transform of Eq. (|17|) gives 

Scxrs(u) = D ac W ca S(uj - gj ca ) 

ac 

+ J2i D ^ W ca}*H^ + ^ca)- (20) 
ac 

As was done in Eq. (fH)]) for RIXS, we can express 
the Doorway and the Window in terms of the Green's 
functions and the CIS coefficients. Details are given in 
the Appendix 151 We then obtain 

D ac := P(a) 1 ^2fii m fi mj f*. ik J ^G$ m) {uj - uj ca ) 
ijk 

x E^Etiu - u ca ) (21) 
W ca := 2^ ///.„//„ ,/,:,, J — lm{G^ n \oj - u ca )} 

ijk 

x E 2 {uo)E*{u J -u Jca ). (22) 

Note that the exact expression for the pump-probe sig- 
nal, Eq. (fTS)) with (|A6|) and (|A7[) , requires the compu- 
tation of many-body states for the TV and the N + 1 
valence electron system. In contrast, the Green's func- 
tion expressions, Eq. (p?Tj) and (J22), only require the 
one-particle Green's functions and valence excited states 
at the CIS level. These may be obtained using standard 
electronic structure computer codes. In Appendix 
we show that Eqs. fA6 |) -([A7 | reduce to (J2TJ) - (J22J) within 
the CIS approximation. 

The physical picture offered by the Doorway and the 
Window in Eq. pi"]) and is shown schematically in 
Fig. 4. The system interacts twice with the pump cre- 
ating a valence excited state electronic wave packet [Eq. 
<|T~8|) together with (f2Tj) ]. This wave-packet evolves freely 
for a time t, which is not limited by the core-hole lifetime. 
The two dipole matrix elements Hi m and fj, m j represent 
the two interactions with the X-ray pulse which create 



and annihilate the core-hole at m. The Green's function 
Gkj(uj) together with coefficients f*c;ik represents the 
relevant dynamics. The Window, Eq. (fT!))) with (J22J), is 
similarly created by the probe X-ray pulse and the signal 
is given by the overlap of these two wavepackets. 



II. CONCLUSIONS 

We have presented a many-body Green's function the- 
ory of stimulated attosecond Raman X-ray scattering. 
The Green's function expression avoids the explicit com- 
putation of the ground and the core-hole excited states of 
the system. The single particle Green's functions can be 
obtained by the self-consistent solution of Hcdin's equa- 
tions [3ll . I32I 32 1 which is implemented in standard com- 
puter codes 19(. We have assumed that the core-hole 
is localized and only enters as a parameter. This al- 



lows to express the signal using single-particle Green's 
functions. The "shakeup" and "shakeoff" excitations 
[19( in which more than one electron is excited due to 
creation of a core-hole is included within the CIS ap- 
proximation. Equation (|12p expresses the correlation 
function of exciton operaotors in terms of the single- 
electron Green's functions and the particle-hole prop- 
agator for valence orbitals. The propagator is subse- 
quently expressed in terms of the CIS expansion coef- 
ficients for the valence excited states. This approxima- 
tion can be relaxed by calculating this propagator using 
the time-dependent Hartree Fock (TDHF). Since TDHF 
is often insufficient, higher-order corrections can be de- 
rived systematically using the many-body expansion for 
the particle-hole propagator [34]. The present formula- 
tion can be extended to incorporate core-hole and nuclear 
dynamics through the corresponding Green's functions 
using eq uation of motion technique [351 ] as was done in 
Ref. M- 
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APPENDIX A: DERIVATION OF EQ. [TT] 

The correlation function expression corresponding to 
the diagrams shown in Figs. 3a and 3b is given by 



/OO ft 
dt / 
-OO J — ! 



dt-, 



dt x {B(t 3 )B\t)B(t 2 )B\t x )) 



dh / dh{B(h)B^(ti)B{h)B^t)) 

> J —CO 

E*(t2 - ^E^h ~ r)E 2 (t)E*(t 3 ). 



(Al) 



Since the pulses are temporally-well-separated and are 
short compared to the time delay r, the upper limits of 
t 2 and t 3 integrals in the first and second terms may be 
safely extended to infinity. Interchanging the integration 
variables t\ and t 2 in the first term we obtain 

Scxrs{t) = 

/OO ft f OO ft\ 

dt I dt 3 / dti / dt 2 E 2 {t)E 2 (t 3 ) 
-oo J — oc J —oo J —oo 

( B{U)B^t)B{h)B\t 2 ) )E*(tx - r)E x (t 2 - t) 

+ (B(t 2 )B\t 1 )B(t 3 )B\t))E* 1 (t 2 -T)E 1 (t 1 -T) . 

(A2) 



where the Doorway wavepacket is given by Eq. (fT8|) with 



D ac = ^P(a)S ae B| c J 



duj Ei{uj)El(uj + LUac), 



2n lu + LU e 



IT] 



(AS) 



The Window wavepacket is similarly given by Eq. 
([T9l)with 



W ca = ^B ce Bl f ^ 

e 

E 2 (lo)E 2 (uj + u ac ). 



1 



1 



lu + uj ec + irj lj + uj ec — irj 

(A7) 



Since in the limit 77 — ^ 

1 1 



lj + lj cc ± irj 



PP- 



LU + LJ e 



T l7T(5(w + LJec), (A8) 



where PP(l/a;) denotes the principal part of l/x, the 
frequency integration in Eq. (|A7|) can be performed, re- 
sulting in 

W ca = 2nJ2B C eBl a E 2 (Lo ce )E* 2 {Lj ac ). (A9) 



APPENDIX B: GREEN'S FUNCTION 
EXPRESSION FOR THE PUMP-PROBE SIGNAL 

We start with Eq. (|A1[) and allow for the creation of a 
core-hole at m and n. The interaction with E\ (El) [E 2 
(E 2 )j creates (destroys) a core-hole at m [n]. We then 
have 

00 ct 



The correlation functions in Eq. (|A2[) can be expanded 
in terms of the many-body states of the system. By 
inserting the identity operator / = ~Y^ V \ v )( v \i where 
v = c, e, we obtain 

(B(t 3 )B\t)B(t 1 )B^(t 2 ))= £ P(a)B ae Bt c B ce ,Bt, a 

acee' 

where B ac = (a\B\c), etc. and P(a) is the equilibrium 
probability of the system to be in state \a). A similar 
expression can be obtained for the other correlation func- 
tion. Substituting these results in Eq. (IA2[) and defining 
the Fourier transform of the pulse envelope as 



/OO ft 
dt / dt 
-OO J — CO 



E(t) 



dw 
2^ 



1 E(lo), 



(A4) 



the signal can be expressed in the Doorway /Window rep- 
resentation as [1, H, [13 ] 



00 J —00 

fOO pt2 

><■ < I dt 2 dt x (B n (t 3 )Bl(t)B m (t 2 )Bl(t x )) 

OO J —CO 

CO pti 

dt x / dt2(B m (t 2 )Bl(h)B n (t 3 )Bl(t)) 

-00 J —CO 

x E{{t 2 -T)E x {t x -T)E 2 (t)E*{t 3 ). (Bl) 

Inserting the projection operators twice in Eq. the 
correlation function in the first term can be expressed as 

(B n (t 3 )Bl(t)B m (t 2 )Bl(t x )) 

« Y, P ( a )( a \ B n {t 3 )V n Bl(t)B m (t 2 )V m Bl{h) \a) 

a 

a ijkl 

x (a\c k clU(t 3 - t)V n c n c]U(t - t^CiC^VmUih - t x )c m c]\a) 

(B2) 



Scxrs{t) = J2i D -cW ca e- 



+ c.c 



(A5) 



Substituting this back in Eq. (|B1[) . and changing vari- 
ables t 2 — ti = s' and t — t 3 = s, the first term can be 
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expressed as 



/oo />oo 
dt I ds 



dU 



ds'e^ t ~ t ^G ( ^\s)G^\s') 



x (x(i'J)\U(t-t 2 )\ X (i,f)) 

x E 2 {t)E*{t-s)E* 1 {t 2 -T)E 1 (t 2 -s'-T) (B3) 
where the \x(hj)) basis is defined in Eq. ((9]) and 

6( S ){a\clc t U( S )c}c m \a) = t6( S )c^ s G[f (a) 

0( S )(a\cl Cl U(~ S )c]c m \a) = -i9(s)e^ s G^ (s\BA) _ 

By inserting tne identity operator, |e)(e|, before and 

Inserting the identity operator in terms of the CIS states after the evolution operator and performing the integra- 
tion over time, Eq. (|C1|) reduces to 



APPENDIX C: EQUIVALENCE OF EQS. 
(fA6 l> - (fA7|) AND ([2l])-([22]) WITHIN THE CIS 
APPROXIMATION 

The pump-probe signal is given by Eq. (|15p where the 
Doorway and the Window are given in terms of the many- 
body states by Eqs. (|A6|) and (|A7|) . Green's function 
expressions are given by Eqs. (|2"Tj) and l|22p. 

We first expand the GF in the many-body states. Us- 
ing Eqs. (JTTJ) and (|B4|) . we can write 

G$ m \u 1 + u, ac ) = 

- i / dte-^-^ialclaU^ctcmla). (CI) 



in (Eq. (T3J)), Eq. (|B3J) finally becomes 

f^imf^mj f^knf^nl fc;li' fc-i 

ijkl i'j' 



ac 

OO (• oo 



x / dt I dse- lu ™ t G t $\s)E 2 {t)E*{t- s) 

-oo Jo 



dt-: 



-co JO 



r(™\ t-^ ( a \cict,\e)(e\c m cl\a) 

G* (m) (u + uiac) = -Y 3 ■ (C2) 

J ^ + LO ec - IT] 

Substituting this in Eqs. (piTj) and (f2"2")l . we obtain 
D M := 



ds'e^^GfflWEifa - r)£i(ta - a' - r) 

(B5) 



£p(a)/i ifn /i mj -/* ifc (a|cfcc^|e)(e|c m ct|a) 
El(w)Ei_(w - uj ca ) 



Performing Fourier transform [Eq. (|A4[) ]. Eq. (IB3 
can also be expressed as 

h = 2Re $>(<*)££ f^imf^mj f^knf^nlfc;li' fc;ij' 
ac ijkl i'j' 

x / - io ca )E 2 (uj)E* 2 (uj - WcB )e*"- T 



W m ■= - 



2n lo + a; ec — irj 
y^Jtkn^ij fc-ji(a\cici\e) (e\c n cj\a) 



(C3) 



ijke 

duo 
2^ 



1 



1 



w + w ec + «?7 w + w ec — irj 

X E 2 {u)E2{u - Uca). 



(C4) 



We next consider the expressions for the Doorway and 
Window in terms of the many-body states. Using Eq. 

as, 



/ ^ G fr ] ^ - Uca)Et(u)Ei(u - uj ca ) (B6) ®> E qs. CM| and (Ml can be expressed 

D ac = ^P(a)fi im fj, mj (a\c i cl n \e)(e\c m c i j \c) 



where Gf^ m \uj) is defined in Eq. (fTTj) . 

Proceeding along the same steps that lead from Eq. 
(|B2|) to (|B6|) . and taking the complex conjugate (this is 
allowed since we are looking for the real part), the second 
term in Equation (|B1[) can be recast as 



duj Ei (oo)E* (lo + io ac ) 
2n lo + lu cc + irj 



(C5) 



and 



h = -2Re^P(a)^^ 

f^imf^mj f^knf^nlfc;lj' /, 

ac ijkl i'j' 

dto 



W ca = Y^kn^nj{c\c k cl L \e){e\c n c\\a) 

1 1 



x / F G 5 n) V - uj ca )E 2 (uj)E*(uj - uj ca )e 



du) 
2^ 



lo + tu ec + irj lo + Lo ec — irj 
E 2 {uj)E;(u + uj ac ). (C6) 

x J -^-G^ m)ji (uj - uj ca )El(uj)Ei(u; - to ca ). (B7) Using the CIS expansion for the excited states |c), Eq. 

(fl"3|) . we can write 



The signal is finally given by 

Scxrs(t) = h + h 
which results in Eqs. (|2ip and 



(B8) 



(e\c m c^\c) = £ f C ;i'f (e\cmc]ci'c], \a) 
i'j' 

= ^h;jj'{e\c m c],\a) 

j' 



(C7) 
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where in going to the second line the sum over i' can be 
done since state |e) corresponds to single excited N + 1- 
valence electrons. Substituting (|C7j) in (|C5|) we get 

£>ac = y^ j P{a)^i m ^ m jfc-jk{a\cic\ n \e){e\c m c\\a) 

ijke 

J 2ir ui + cj ec + i?7 



Taking the complex conjugate and interchanging dummy 
indices i and j, Eq. (|C8p becomes same as (|C3p . Simi- 
larly one can show that Eq. (|C6[) is same as (|C4|1 within 
CIS approximation. 

This proves the equivalence of Eqs. (|A6[) - (|A7[) and 
Eqs. (j2T|)-(i22l). 



[1] M. Drescher, M. Hentschel, R. Kienbereger, G. Tempea, 
C. Spielman, G. A. Reider, P. B. Corkum and F. Krausz, 
Science 291, 1923 (2001). 

[2] R. F. Service, Science 298, 1356 (2002). 

[3] Attosecond Spectroscopy: Special Issue of Science 317, 
August (2007). 

[4] I. Thomann, E. Gregonis X. Liu R. Trebino, A. S. 

Sandhu, M. M. Murnane and H. C. Kapteyn, Phys. Rev. 

A 78, 011806(R) (2008); L. Miaja-Avila, G. Saathoff, S. 

Mathias, J. Yin, C. La-o-vorakiat, M. Bauer, M. Aeschli- 

mann, M. M. Murnane and H. C. Kapteyn, Phys. Rev. 

Lett. 101, 046101 (2008). 
[5] F. Dorchies, F. Blasco, C. Bonte, T. Caillaud, C. Four- 

ment and O. Peyrusse, Phys. Rev. Lett. 100, 205002 

(2008). 

[6] S. Mukamel, Phys. Rev. B 72, 235110 (2005) 

[7] V. Felicissimo, F. Guimaraes, F. Gel'mukhanov, A. Cesar 

and H. Agren, J. Chem. Phys. 122, 094319 (2005). 
[8] I. V. Schweigert and S. Mukamel, Phys. Rev. A 76, 

012504 (2007). 

[9] S. Tanaka and S. Mukamel, Phys. Rev. Lett. 89, 043001 
(2002). 

[10] C. Bressler and M. Chergui, Chem. Rev. 104, 1781 
(2004). 

[11] D. A. Oulianov, I. V. Tomov, A. S. Dvornikov and P. M. 
Rentzepis, PNAS, 99, 12556 (2002). 

[12] L. X. Chen, X. Zhang, E. C. Wasinger, K. Attenkofer, 
G. Jennings, A. Z. Muresan and J. S. Lindsey, J. Am. 
Chem. Soc. 129, 9616 (2007); L. X. Chen, Ann. Rev. 
Phys. Chem. 56, 221 (2005). 

[13] R. Santra, C. Buth, E. R. Peterson, R. W. Dunford, E. P. 
Kanter, B. Krassig, S. H. Southworth, L. Young, Journal 
of Physics: Conference Series 88, 012052 (2007). 

[14] S. Mukamel, Principles of Nonlinear Optical Spec- 
troscopy, Oxford University Press, New York (1995). 

[15] C. A. Ashley and S. Doniach, Phys. Rev. B 11, 1279 
(1975). 

[16] J. J. Rehr and A. L. Ankudinov, Coordination Chem. 
Rev. 249, 131 (2005). 

[17] J. A. Soininen, A. L. Ankudinov and J. J. Rehr, Phys. 
Rev. B 72, 045136 (2005). 

[18] A. L. Ankudinov, B. Ravel, J. J. Rehr and S. D. Conrad- 
son, Phys. Rev. B 58, 7565 (1998). 

[19] J.J. Rehr and R. C. Albers, Rev. Mod. Phys. 72, 621 
(2000). 



[20] Yejun Feng, J. A. Soininen, A. L. Ankudinov, J. O. Cross, 

G. T. Seidler, A. T. Macrander, J. J. Rehr, and E. L. 
Shirley, Phys. Rev. B 77, 165202 (2008). 

[21] F. Gel'mukhanov and H. Agren, Phys. Rep. 312, 87 
(1999). 

[22] M. Magnuson, L.-C. Duda, S. M. Butorin, P. Kuiper and 
J. Nordgren, Phys. Rev. B 74, 172409 (2006); C. Sathe, 
F. F. Guimaraes, J.-E. Rubensson, J. Nordgren, A. Agui, 
J. Guo, U. Ekstrom, P. Norman, F. Gel'mukhanov, and 

H. Agren, Phys. Rev. A 74, 062512 (2006). 

[23] F. de Groot and A. Kotani, "Core Level Spectroscopy of 
Solids", CRC Press (Taylor & Francis Group), New York 
(2008). A. Kotani and S. Shin, Rev. Mod. Phys. 73, 203 
(2001). 

[24] H. Sternemann, J. A. Soininen, C. Sternemann, K. 
Hamalainen and M. Tolan, Phys. Rev. B 75, 075118 
(2007); S. Galambosi, M. Knaapila, J. A. Soininen, 
K. Nygard, S. Huotari, F. Galbrecht, U. Scherf, A. 
P. Monkman and K. Hamalaineen, Macromolecules 39, 
9261 (2006). 

[25] Y. Hikosaka, Y. Velkov, E. Shigemasa, T. Kaneyasu, Y. 

Tamenori, J. Liu and F. Gel'mukhanov, Phys. Rev. Lett. 

101, 073001 (2008). 
[26] L. E. Fried and S. Mukamel, J. Chem. Phys. 93, 3063- 

3071 (1990). 

[27] C. Marx, U. Harbola and S. Mukamel, Phys. Rev. A 7, 
022110 (2008). 

[28] A. Fetter and J. D. Walecka, Quantum Theory of Many- 
Particle Systems, Dover Publication, New York (2003). 

[29] H. Haug, A. -P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors, Springer- Verlag, Berlin, 
Heidelberg, (1996) 

[30] T. Privalov, F. Gel'mukhanov and H. A gren, Phys. Rev. 
B 64,165116 (2001). 

[31] L. Hedin, Phys. Rev. 139, 796 (1965). 

[32] F. Aryasetiawan and O. Gunnarsson, Rep. Prog. Phys. 
61, 237 (1998). 

[33] G. Onida, L. Reining and A. Rubio, Rev. Mod. Phys. 74, 
601 (2002). 

[34] J. Lindenberg and Y. Ohrn, Propagators in Quantum 
Chemistry, 2nd ed.,John Willy & Sons, Inc., New Jersy 
(2004). 

[35] D. C. Lengreth, Phys. Rev. 182, 973 (1969). 



